WORLD INTELLECTUAL PROPERTY ORGANIZATION 
International Bureau 




PCX 

INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(51) International Patent Classification ^ : 

H04B 7/OOS, H04L 25/03, H04B 1/707 



A2 



(11) International PubUcation Number: WO 98/44655 

(43) International Publication Date: 8 October 1998 (08.10.98) 



(21) International Application Number: PCT/SE98/00518 

(22) International Filing Date: 20 March 1998 (20.03.98) 



(30) Priority Data: 

08/827,169 



27 March 1997 (27.03.97) 



US 



(71) Applicant: TELEFONAKTIEBOLAGET LM ERICSSON 

(publ) [SE/SE]; S-126 25 Stockholm (SE). 

(72) Inventors: ANDERSSON, SOren; Bergsatersvagen 9D, S-191 

40 Sollentuna (SE). PELIN. Per; Odinsgatan 22B, S-41 1 

03 Gttteborg (SE). 

(74) Agent: ERICSSON RADIO SYSTEMS AB; Common Patent 
Dept., S-164 80 Stockholm (SE). 



(81) Designated States: AL, AM, AT, AU, AZ, BA, BB, BG, BR, 
BY, CA, CH, CN, CU, CZ, DE. DK. EE. ES. FI, GB, GE, 
GH. GM. GW, HU. ID, IL, IS, JP. KE, KG, KP, KR, KZ, 
LC, LK, LR, LS, LT. LU, LV, MD, MG, MK, MN, MW, 
MX, NO, NZ, PL, PT, RO, RU, SD, SE, SG, SI. SK, SL, TJ, 
TM, TR, TT, UA, UG, UZ, VN, YU, ZW, ARIPO patent 
(GH, GM, KE, LS. MW, SD, SZ, UG. ZW), Eurasian patent 
(AM, AZ, BY, KG, KZ, MD, RU , TJ, TM), European patent 
(AT, BE, CH, DE, DK, ES, FI, FR, GB, GR, IE; IT, LU, 
MC, NL, PT, SE), OAPI patent (BF, BJ. CF, CG, CI, CM, 
GA, GN, ML, MR, NE, SN, TD, TG). 



Published 

Without international search report and to be republished 
upon receipt of that report. 



(54) Title: SELF-SYNCHRONIZING EQUALIZATION TECHNIQUES AND SYSTEMS 

/- 

(57) Abstract 

Methods and systems for processing received radio signals are described wherein 
the finite alphabet quality of digital transmissions is utilized to improve performance. 
Nonsynchronous sampling of received signals introduces intersymbol interference which 
is compensated for by these methods and systems. Different types of signal modulation, 
and their impact on this type of intersymbol interference are discussed. A robust diversity 
combining technique, usable in conjunction with antenna arrays is developed. 







TRA// 
SBQUi 
TO OEM 

r» 


40WN 
//A/0 
=NCE 
ERATE 



USer„(t)ANO 
LS-METHOD TO 
CALCULATE 
INITIAL 
CHANNEL ESVMATE 






BNSOOCID: <WO_fl844655A2j_> 




FOR THE PURPOSES OF INFORMATION ONLY 
Codes used to identify States party to the PCT on the front pages of pamphlets publishing international applications under rtie PCT. 



AL 


Albania 


ES 


Spain 


LS 


Lesotho 


SI 


Stovenia 


AM 


Armenia 


FI 


Finland 


LT 


Lithuania 


SK 


Slovakia 


AT 


Austria 


FR 


France 


LU 


Luxembourg 


SN 


Senegal 


AU 


Australia 


GA 


Gabon 


LV 


Latvia 


sz 


Swaziland 


AZ 


Azerbaijan 


GB 


United Kingdom 


MC 


Monaco 


TD 


Chad 


BA 


Bosnia and Herzegovina 


GE 


Georgia 


MD 


Republic of Moldova 


TG 


Togo 


BB 


Barbados 


GH 


Ghana 


MG 


Madagascar 


TJ 


Tajikistan 


BE 


Belgium 


GN 


Guinea 


MK 


The former Yugoslav 


TM 


'nirkmenistan 


BF 


Burkina Paso 


GR 


Greece 




Republic of Macedonia 


TR 


Turkey 


BG 


Bulgaria 


HU 


Hungaiy 


ML 


Mali 


TT 


Trinidad and Tobago 


BJ 


Benin 


IE 


Ireland 


MN 


Mongolia 


UA 


Ukraine 


BR 


Brazil 


IL 


Israel 


MR 


Mauritania 


UG 


Uganda 


BY 


Belarus 


IS 


Iceland 


MW 


Malawi 


US 


United States of America 


CA 


Canada 


IT 


Italy 


MX 


Mexico 


UZ 


Uzbekistan 


CF 


Central A^can Republic 


JP 


Japan 


NE 


Niger 


VN 


Viet Nam 


CG 


Congo 


KE 


Kenya 


NL 


Netherlands 


YU 


Yugoslavia 


CH 


Switzerland 


KG 


Kyrgyzstan 


NO 


Norway 


ZW 


Zimbabwe 


CI 


C6tc d'lvoire 


KP 


Democratic People's 


NZ 


New Zealand 






CM 


Cameroon 




Republic of Korea 


PL 


Poland 






CN 


China 


KR 


Republic of Korea 


PT 


Portugal 






CU 


Cuba 


KZ 


Kazakstan 


RO 


Romania 






CZ 


Czech Republic 


LC 


Saint Lucia 


RU 


Russian Federation 






DE 


Germany 


LI 


Liechtenstein 


SD 


Sudan 






DK 


Denmaifc 


LK 


Sri Lanka 


SE 


Sweden 






EE 


Estonia 


LR 


Liberia 


SG 


Singapore 







BNSDOCID: <WO_9844655A2J_> 



wo 98/44655 PCT/SE98/00518 



SELF-SYNCHRONIZING EQUALIZATION TECHNIQUES AND SYSTEMS 

BACKGROUND 

In recent years, wireless communication systems have been used to 
convey a variety of information between multiple locations. With digital 
5 communications, information is translated into a digital or binary form, .referred 
to as bits, for communications purposes. The transmitter maps this bit stream 
into a modulated symbol stream, which is detected at the digital receiver and 
mapped back into bits and information. 

In digital wireless communications, the radio environment presents many 

10 difficulties that impede successful communications. One difficulty is that the 
signal level can fade, because the signal may travel in multiple paths. As a 
result, signal images arrive at the receiver antenna out of phase. This type of 
fading is commonly referred to as Rayleigh fading or fast fading. When the 
signal fades, the signal-to-noise ratio becomes lower, causing a degradation in 

15 the cormnunication link quality. 

A second problem occurs when the multiple signal paths are much 
different in length. In this case, time dispersion occurs in which multiple fading 
signal images arrive at the receiver antenna at different times, giving rise to 
signal echoes. This causes intersymbol interference (ISI), where the echoes of 

20 one symbol interfere with subsequent symbols. 

Raleigh fading can be mitigated by using diversity, such as antetina 
diversity, at the receiver. The signal is received on a plurality of antennas. 
Because the antennas have slightly different locations and/or antenna patterns, the 
fading levels on the antennas are different. In the receiver, these multiple 

25 antenna signals are combined either before or after signal detection using such 
techniques as maximal-ratio-combining, equal-gain-combining, and selective 
combining. These techniques are well known to those skilled in the an and can 
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be found in standard textbooks, such as W.C.Y. Lee, Mobile Communications 
Engineering, New York: McGraw-Hill, 1982. 

The tune dispersion can be mitigated by using an equalizer. Common 
forms of equalization are provided by linear equalizers, decision-feedback 
5 equalizers, and maximum-likelihood sequence-estimation (MLSE) equalizers. A 
linear equalizer tries to undo the effects of the channel by filtering the received 
signal. A decision-feedback equalizer exploits previous symbol detections to 
cancel out the intersymbol interference from echoes of these previous symbols. 
Finally, an MLSE equalizer hypodiesizes various transmined symbol sequences 
10 and, widi a model of die dispersive channel, determines which hypodiesis best 
fits the received data. These equalization techniques are well known to those 
skilled in die art, and can be found in standard textbooks such as J.G. Proakis, 
Digital Communications, 2nd ed. New York: McGraw-Hill, 1989. 

Of the three common equalization techniques, MLSE equalization has 
15 been considered preferable from a performance point of view. In die MLSE 
equalizer, all possible transmitted symbol sequences are considered. For each 
hypothetical sequence, die received signal samples are predicted using a model of 
die multipadi chaimel. The difference between die predicted received signal 
samples and die actual received signal samples, referred to as die prediction 
20 error, gives an indication of how good a particular hypodiesis is. The squared 
magnitude of die prediction error is used as a metric to evaluate a particular 
hypothesis. This metric is accumulated for different hypodieses for lise in 
determining which hypodieses are better. This process is efficiendy realized 
using die Viterbi algoridim, which is a form of dynamic programming. 
25 Ideally, die diversity combining process and the equalization process 

should be combined in some optimal way. Recent research has shown diat for 
MLSE equalization, diversity combining should be done widiin the equalizer. 
This research can be found in W.H. Sheen and G.L. Stflber, "MLSE equalization 
and decoding for muitipadi-fading channels," IEEE Trans. Commun., vol. 39. 
30 . pp. 1455-1464. Oct. 1991; Q. Liu and Y. Wan "An adaptive maximum- 
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likelihood sequence estimation receiver with dual diversity combining/selection, " 
Ind. Symp. on Personal, Indoor and Mobile Radio Commun,, Boston, Mass,, pp, 
245-249, Oct. 19-21, 1992, and Q. Liu and Y. Wan, "A unified MLSE detection 
technique for TDMA digital cellular radio," 43rd IEEE Vehicular Technology 
Conference, Seacaucus, NJ., pp. 265-268, May 18-20, 1993. In the above 
mentioned research, diversity combining is performed by adding together the 
magnitude squared prediction errors from different diversity channels when 
forming metrics. 

The use of antenna arrays at base stations in a mobile communication 
systems has also been proposed as a technique for increasing capacity and 
performance. The most common approach for processing the information 
gatiiered by each antenna associated with a particular signal is based on direction 
of arrival (DOA) estimation followed by beamforming, i.e. combining the vector 
signal from the array to a scalar signal (spatial filtering), before detection. 
However, this approach does not fully exploit the spatial stmcture of the channel. 
A better way is to use an algorithm that is adaptive in the spatial domain and 
which also takes the quality that the transmitted signal has a finite alphabet (e.g., 
O's and I's) into account. Examples of such algorithms are the recently 
proposed iterative least squares with projections (ILSP) algorithm and the 
decoupled weighted least squares with projections (DWILSP) algorithm. The 
decoupled algorithm is similar to ILSP in performance, but is computationally 
cheaper. 

Both ILSP and DWILSP are, in their original formulation, limited to use 
on frequency-flat (i.e., non time-dispersive) channels. However, in many mobile 
conmiunication systems, the channel cannot be modelled as frequency-flat. To 
treat time-dispersive channels, extensions to the iterative least squares approaches 
have also been presented. These algorithms are unformnately quite complex, 
both regarding computational aspects and detection procedures involved. 

Another drawback associated with these conventional algorithms is their 
requirement of precise synchronization. Although the DWILSP algorithm can be 
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used to process signals received from unsynchronized cochannel users, 
synchronization with the signal of interest is stUl assumed, i.e.. the signal of 
interest is assumed to be sampled correctly in accordance with the symbol 
timmg. In practice, this assumption may not hold true, since perfect symbol 
5 timing is difficult to achieve. For example, in certain types of systems, e.g., 
time division multiple access (TDMA) systems which use short transmission 
bursts, proper sample timing is extremely difficult to guarantee. Thus, as wiU be 
illustrated in the simulations performed by Applicants and described below, the 
conventional DWILSP algorithm suffers significant degradation (e.g., increased 
10 bit error rate) when timing errors are introduced mto the sampled signal. 

Several techniques have been proposed which use oversampling, i.e., 
taking more than one time discrete sample during each symbol interval, to handle 
the problems associated with unsynchronized signals. The DWILSP algorithm, 
however, is designed to use only one sample per symbol interval and, therefore, 
15 is not amenable to these types of solutions. 

Accordingly, it would be desirable to provide a technique for estimating 
symbols using the DWILSP algorithm fi-om unsynchronized signals sampled at 
the symbol rate. Moreover, it would also be desirable to use the DWILSP 
algorithm to obtain improved diversity combining. 

20 SUMMARY 

According to exemplary embodiments of the present invention, these and 
other drawbacks and problems associated with the conventional DWILSP 
algorithm, and similar techniques for processing received radio signals, are 
overcome by providing self-synchronizing techniques which provide improved 

25 performance for nonsynchronously sampled signals. For example. Applicants 
have recognized that nonsynchronously sampled signals create additional 
intersymbol interference (ISI) which should be compensated for in order to 
improve detection performance. This additional ISI is different than that 
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described above in that it is parameterizable (and therefore readily determinable) 
based upon timing error and modulation type. 

Thus, exemplary embodiments of the present invention teach the provision 
of compensation schemes which, for example, modify the conventional DWILSP 
technique to compensate for the ISI introduced by nonsynchronous sampling. A 
specific example is given for MSK modulation, although the present invention 
can be applied to any type of modulation with adaptations which will be apparent 
to those skilled in the art. 

In addition to providing self-synchronizing processing techniques, 
exemplary embodiments of the present invention also provide for robust diversity 
combining which outperforms conventional techniques, e.g., RAKE diversity 
combining. By using the DWILSP technique to provide temporal combining of 
spatio-temporal signal estimates created using an adapted version of the RAKE 
algorithm, exemplary embodiments of the present invention are able to 
significantly improve upon prior diversity combining techniques. 



The feamres and advantages of Applicants' invention will be understood 
by reading this description in conjunction with the drawings, in which: 

FIG. 1 is a block diagram of an exemplary cellular radio telephone 
20 system in which the present invention may be applied; 

FIG. 2 illustrates an exemplary antenna array and processing strucmres 
associated therewith; 

FIG. 3 is a flowchan depicting an exemplary self-synchronizing technique 
according to the present invention; 
25 FIG. 4 is a graph illustrating simulation results in terms of bit error rate 

for BPSK modulated signals processed according to both the conventional 
DWILSP technique and self-synchronizing techniques according to the present 
invention; 



BRIEF DESCRIPTION OF THE DRAWINGS 
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FIG. 5 is a graph illustrating simulation results in terms of root mean 
square delay for BPSK modulated signals processed according to self- 
synchronizing techniques according to the present invention; 

FIG. 6 is a graph illustrating simulation results for MSK modulated 
signals processed according to both the conventional DWILSP technique and a 
self-synchronizing technique according to the present invention; 

FIG. 7 is a graph illustrating simulation results for GMSK modulated 
signals processed according to both the conventional DWILSP technique and a 
self-synchronizing technique according to the present invention; 

FIG. 8 is a block diagram of a conventional RAKE combiner; 

FIG. 9 is a block diagram of another known RAKE combiner using the 
DWILSP technique to provide signal estimates; 

FIG. 10 is a block diagram of a RAKE combiner according to an 
exemplary embodiment of the present invention; 

FIG. 11 is a block diagram of a RAKE combiner according to an another 
exemplary embodiment of the present invention; 

FIG. 12 is a flowchart illustrating steps associated with an exemplary 
diversity combining technique according to the present invention; 

FIG. 13 is a graph illustrating the results of a first simulation used to 
demonstrate the performance of an exemplary diversity combining technique 
associated with the present invention; and 

FIG. 14 is a graph illustrating the results of a second simulation used to 
demonstrate the performance of an exemplary diversity combining technique 
associated with the present invention. 

DETAILED DESCRIPTION 
The following description is scripted in terms of a cellular 
radiocommunication system, but it will be understood that Applicants* invention 
is not limited to that environment. To provide an understanding of various 
exemplary receivers and systems within which structures and techniques 
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according to the present invention can be implemented, the following summarizes 
an exemplary cellular radiocommunication system. 

FIG. 1 is a block diagram of an exemplary cellular radiocommunication 
system, including an exemplary base station 110 and mobile station 120. The 
5 base station includes a control and processing unit 130 which is connected to the 
mobile switching center (MSG) 140 which in turn is connected to the PSTN (not 
shown). General aspects of such cellular radiocommunication systems are known 
in the an, as described by the above-cited U.S. patent applications and by U.S. 
Patent No. 5,175,867 to Wejke et aL, entitled "Neighbor- Assisted Handoff in a 

10 Cellular Communication System," and U.S. Patent Application No. 07/967,027 
entitled "Multi-Mode Signal Processing," which was filed on October 27, 1992, 
both of which are incorporated in this application by reference. 

The base station 1 10 handles a plurality of traffic channels through a 
traffic channel transceiver 150, which is controlled by the control and processing 

15 unit 130. Also, each base station includes a control channel transceiver 160, 

which may be capable of handling more than one control channel. The control 
channel transceiver 160 is controlled by the control and processing unit 130, 
The control channel transceiver 160 broadcasts control information over the 
control channel of the base station or cell to mobiles locked to that control 

20 channel. It will be understood that the transceivers 150 and 160 can be 

implemented as a single device, like the traffic and control transceiver 170 in the 
mobile station, for use with control channels and traffic channels that share the 
same radio carrier frequency. The traffic channels can be used in a dedicated, 
connection-oriented maimer to transmit information, e.g., for a voice connection, 

25 where each channel is used continuously for a period of time to support 

transmission of a single stream of information or in a packet-oriented manner 
where each channel can be used to send independent units of information 
associated with different information streams. 

Transceivers 150 and 160 may have dedicated antennas 170 and 180 

30 which, using a duplex filter, transmit and receive signals for processing therein. 
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Altematively, base station 110 may be provided with an antenna array as 
depicted in FIG. 2. The antenna array will have some rnmiber m of antenna 
elements 200, where m> =2. Each signal creates a response on each antenna 
element 200, which response is processed (e.g., filtered, downconverted, etc) in 
5 receive processing blocks 210. The processed signal responses are used to 

generate a channel estimate and a signal estimate Sk(t) for each sampling time 
instance i as shown in blocks 220. The manner in which these estimates are 
created and combined are described below with respect to exemplary 
embodiments of the present invention. 

^® In order to have a complete understanding of the present invention, it is 

first beneficial to consider the origins thereof, in particular the DWILSP 
technique referred to above. A similar description of this conventional technique 
can be found in die article entitled "Decoupled Separation of Digitally Modulated 
Signals Arriving at an Antenna Array", authored by P. Pelin et al., published in 

15 Proc. of RVK 96, Lulea, Sweden, June 1996, the disclosure of which is 
expressly incorporated here by reference. 

The Conventional DWTLSP Technique 

In an environment with multipath propagation, the output of an m-element 
array can be expressed as: 

d <lk 

^(^) = EE«(9ju)ViuP;^i('-f*/) * v(f) (1) 
i=i t-i 

where d is the number of signals impinging on the array, j^. is the signal from the 
20 k:th user (widi symbols belonging to a finite alphabet) and and is the 

attenuation and time-delay for each of the q„ subpaths. 

Herein, a narrow-band assumption is imposed (i.e., the propagation 
25 delays associated with multipath are much smaller than the inverse bandwidth of 

the signals), so that 5,^(1-7^) expi-JwrMO - Equation (1) can tiius be rewritten 

as: 
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d (2) 

^(0 = Y^^iPi^k^t) + v(t) 



where ^ f=i yjcP^vH^^nX^k) , called the spatial signature, is the sum 
of multipath array responses due to signal k. Assuming the d signals are symbol- 
synchronized, which assumption (as mentioned above) increases the bit error rate 
5 associated with the received . signal under many practical conditions, the antenna 
outputs are passed through a filter matched to the transmit pulse, and sampled at 
the symbol rate /J = 1/T to yield the corresponding discrete-time model: 

x(n) = As(n) + v(n) (3) 

where A^,^ is the collection of total array response vectors, scaled by the signal 
10 amplitudes, i.e., A = [p^a^...p/ij, s(n) = [b^in)...b/n)f, bfn) = ±1, and v(/z) 
is spatially and temporally white noise, A block formulation is obtained by 
taking N snapshots, yielding: 

X(N) = \S(N) + W(N) (4) 
where X„,;/N) = [x(l)..,x(N)]. S^,;/N) = [s(l)...s(N)], and V^,^N) = 
15 [v(l),..v(N)]. The spatial structure of the data is represented by A, while the 
matrix S represents the temporal structure. The above formulation is valid for 
BPSK (binary phase shift keying) signals, but extension to arbitrary linear 
modulation schemes is straightforward. 

By defining one signal (at a time) to be signal of interest (SOI), equation 
20 (4) can be rewritten in the following way: 

XOT-'-AE"'''*^'^ (5) 

where the first signal is taken to be the SOI, without loss of generality. The 
term J(N) thus corresponds to interfering signals plus noise. Since it is desired 
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to estimate the signals with little or no spatial knowledge, a and s can be 
iteratively estimated, based on the formulation in equation (5). 

Given an initial estimate of a spatial signature a, the following weighted 
least-squares criterion function is iteratively minimized: 

^(X-as)'PF(X-as)=^| |PF^(X^as)| 1^ (6) 

5 Here, W should ideally be chosen as y , which can be interpreted as a 

pre whitening of the data vector x(n). However, it can be shown using the matrix 
inversion lemma, that using the inverse of the sample estimate of the covariance 
of the array output produces asymptotically equivalent signal estimates. Equation 
10 (6) can thus be reformulated as follows: 

"^WZ-bsW' (7) 

with Z = ^'^Xand b = ^^o.- For fixed b, the solution to equation (7) with 
respect to s is : 

s = (Jb'by'b'Z = —^S'Z (8) 

WW 

15 Exploiting the finite-alphabet property, this solution is projected onto its closest 
discrete values in the signal space. In the case of BPSK signals, this. projection 
is equivalent to taking the sign of each component in s. The (modified) spatial 
signature b is then updated by minimizing equation (7) with respect to b. The 
solution is: 

b = Zs-(ii V = — 

N (9) 



20 



Note that equation (9) is a temporally matched filter to the current signal 
estimate, whereas (8) represents a spatially matched filter. The process is 
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repeated until s converges, after which the algorithm continues with the next 
signal. 

As mentioned above, the conventional DWILSP algorithm does not take 
into consideration that the symbol sampling is imperfect. Accordingly, the 
5 present invention modifies the aforedescribed technique to handle mtersymbol 

interference caused either by non bit-synchronized sampling or by the modulation 
technique used to process the original signal for transmission over the air 
interface. These modified techniques accordmg to the present invention are 
referred to herein as "self-synchronized" techniques. 

10 Self-Svnchronized Techniques 

Sampling a signal in an unsynchronized manner means, for most 
modulation formats, that intersymbol interference (ISI) is introduced. This fomi 
of ISI is quite different from the ISI caused by a time dispersive propagation 
channel. The reason for this is that ISI caused by unsynchronized sampling has 

15 an underlying structure, i.e., the ISI can be parameterized by the timing error. 
The parameterization of this strucmred kind of ISI differs between 
modulation formats. Therefore, the m9difications made to the DWILSP 
technique according to exemplary embodiments of the present invention will also 
depend on the modulation format. 

20 The effects of ISI due to nonsynchronous sampling are reflected in the 

data model by a modification of the source signal description as: 



d 



(10) 



Here, the ISI is parameterized in the scalar signal Sis,.i(n), and the 
characterization of this ISI depends on the modulation format. In some cases, 
25 there is no ISI at all, for example MPSK modulation with a rectangular pulse 
shape, sampled directly at the symbol rate without a preceding matched filter. 
Nevertheless, in most cases, sampling nonsynchronously leads to ISI, as for 
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example when a signal modulated by minimum shift keying (MSK) is 
nonsynchronously sampled. 

An MSK signal is most often received by direct sampling at the symbol 
rate, without any matched filter, as in the European GSM system and systems 
operating in accordance with the GSM standard. Generally, the received signal, 
nonsynchronously sampled, can be expressed as: 



S(t-1) • COS 



ft 

2T 



(11) 



where T, is the sampling interval and r e [0.1] is the timing error in the sampling 
(relative to T,). From this equation it can be seen diat for an MSK-signal, the 
signal s,si(t) is characterized by: 

a constant envelope, that is, it has the same power mdependent of 
the value for r; 

• the ISI-components are separated in quadramre, where: 
•^(') • ^"1^^] in-phase component, and 



-s{t-l) ■ 



TZ 

2T 



is the quadrature component. 



These qualities make it possible to write the output. x(t). from the m-element 
array antenna for a flat fading channel or a spatio-temporal diversity path as: 



x(t) = h ■ [cti-J-a^] 



sit) 



+ J 



(12) 
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10 



15 



where. 



(13) 



aj = cos 



(14) 



x(t) € C""', and h e C""' describes the propagation channel for a flat fading 
channel (i.e., without time dispersion) or a spatio-temporal diversity path. 
Finally, J denotes any modeling error. To provide even better performance, a 
pre-whitening process can be applied. Pre-whitening is achieved by computing 
the following new quantities. First, the estimated array covariance matrix is 
defined by (with "H" denoting the Hermitean transpose operator): 

= f:x(f)x^(f) ^^^^ 

and then the pre- whitened array output data and channel response vector as, 
resp.. 



and. 



(16) 



(17) 



The self-synchronizing technique according to the present invention for 
detecting/estimating ISI in a single diversity path can now be outlined as follows. 
The flowchart of FIG. 3 provides a visual guide to the below described steps 
according to the present invention. 
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Assume, at block 300, an initial timing error, e.g., t=0.5, and create a 
corresponding signal r,s,(t), using a known training/reference sequence r(t) (which 
is contained as a part of the original finite alphabet signal, s(t), transmitted from 
a mobUe station). Those skilled in the art will appreciate that different systems 
5 provide different known reference sequences in their transmission bursts. For 
example, the GSM system provides a training sequence having 26 bits. 

To continue for the general case, generate, for t=t, up to t=t2, the 
construction denoted r,si(t): 

^/s/O = rit) - j • r(t-l) (18) 

10 the length of this construction will depend on the actual length of the panicular 
training sequence considered. Use this construction, together with the well- 
known Least-Squares (LS) method for parameter estimation, to find an initial 
estimate, g, of the channel response vector at block 310 using the below data 
model (with t=t,...t2): 



15 



z(0 = g ' r,,j(t) (19) 



Having found an initial estimate, the process continues iteratively beginning with 
an estimation of the sampled ISI signal, s,s,(t), employing the LS-method using 
the received pre-whitened data, z(t), and the estimated channel response vector, 
g. as indicated at block 320. Next, the model for the estimated received data can 
20 be rewritten as: 









-//7i(i,^r+l)) 




"2 



'^'^ (20) 



where s(l) is the originally transmitted finite alphabet signal by a mobile station. 
In equation (20). die variables a„ a. and s(t) can then be solved for using the 
conventional DWILSP technique. The relative sampling instance, r, can be 
25 estimated from a, and 
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The next step is to compute an updated estimate of the sampled ISI signal 
according to: Sj^j = a ^ -s(t) -J-a^-sit- 1) 



Then an updated channel response vector, g, can be computed, block 330, 
using the LS-method on the data model: 



^(0 = S ' Sjs^t) 



(21) 



10 



15 



20 



If s(t) has converged as determined at block 340, then the process can be 
terminated, otherwise another iteration begins at step 320. 

Other linear modulation formats, e.g., BPSK, lead to models similar to 
the one presented above for MSK. Note, however, that the ISI parameter vector 
may or may not be a linear function of timing error T; for these other models. 

To test the performance of processing techniques according to the present 
invention, a simulation was conducted that compares the present invention with 
the conventional DWILSP algorithm for signals using BPSK or Gaussian MSK 
modulation. The test simulated a 5-element antenna array that receives two 
signals from nominal DOA:s of [-15°.20°]. The signals are transmitted in bursts 
corresponding to the normal GSM burst, i.e., 148 bits, including a 26 bit 
training sequence in the central part, and three known tail bits at each end. The 
channel was modelled as flat-fading and the scattering cluster width a was 3". 
To simulate Rayleigh fading, independent channel vectors were used for each 
transmitted burst. The average Eb/N, at each antenna-element was set to 5dB. 

In the BPSK case, the performance of the original DWILSP algorithm 
was compared to the self synchronizing technique according to the present 
invention. In the simulation, the self-synchronizing technique was tested twice, 
once using the LS-approach, and a second time using Viterbi equalization to 
facilitate a performance comparison. The timing error introduced by 
nonsynchronous sampling was varied, giving the results shown in FIG. 4. In this 
figure, bit error rate is plotted against timing error. Throughout these 
simulations the foUowme conventions are used. The dashed line reoresents the 
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results for the conventional DWILSP technique, the results for the self- 
synchronizing technique (LS-approach) is shown as a dotted line and the results 
for the self-synchronizing technique (Viterbi approach) is shown using a solid 
line. 

In FIG. 4, it can be seen that either implementation of the present 
invention provides improved perfonnance as compared with the conventional 
DWILSP technique due to its assumption of synchronized bit sampling. The 
numerical problems involved in the LS implementation for timing errors in the 
vicinity of r=0.5 can be seen by the spike in the BER. Using the Viterbi 
algorithm also leads to a performance degradation for tj^O and r?i 1, but this is 
a consequence of the signal power loss involved, and not the Viterbi algorithm 
itself. In some signal processmg applications, for example radar and positioning, 
the timing error t is of more importance than the BER. FIG. 5 shows the root 
mean square (RMS) error of the delay estimate for the LS and Viterbi 
15 implementations of the present invention. 

Repeating the simulation described above but with MSK and GMSK 
modulation for the transmitted signals provides the results shown in FIGS. 6 and 
7, respectively. For both types of modulation, the present invention again 
outperforms the DWILSP technique. For the GMSK case, it can be seen that the 
perfonnance of the self-synchronizing technique according to the present 
invention is only slightly dependent on the timing. The best perfonnance is 
obtained by sampling the received signal between symbol transitions, i.e. t=0.5, 
since the GMSK waveform is much closer to MSK at these mstants. 

As can be seen from the foregoing, the conventional DWILSP algorithm 
required the signal of interest to be sampled con-ectly or the BER perfonnance 
wUl be degraded. By way of contrast, exemplary techniques accordmg to the 
present mvention provides improved perfonnance across the spectrum of timing 
errors and, accordmgly, pennit die signal of interest to be sampled 
nonsynchronously. For some modulation fonnats. some perfonnance degradation 
is introduced, whereas for others, there is no perfonnance degradation involved. 



20 



25 



30 
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The self-synchronizing techniques according to the present invention also 
provide an estimate of the tuning error, either explicitly, or as a function value 
thereof. For example. Equation (13) can be rewritten to provide an estimate of 
the timing error t as the following function value: 



where is here a known quantity, and d, has been estimated by the 
conventional DWILSP algorithm. As a consequence, the self-synchronizing 
version of the DWILSP algorithm can be used for other applications than 
communications, for example radar and positioning. 

Diversity Combining 

The foregoing exemplary embodiments dealt with ISI caused by 
modulation and/or unsynchronized sampling. This type of ISI is deterministic 
since there is a strict underlying parameterization. The following exemplary 
embodiments relate to ISI caused by the propagation channel which, as described 
above, is quite different in namre. In the same way that the mobile radio 
channel spreads the transmitted energy in the spatial domain, i.e., in a stochastic 
manner, the time-dispersion of the channel also causes a spreading of energy in 
the temporal domain. 

Cancelling the effect of the channel dispersion is, as described above, a 
classical problem known as equalization. Conventional techniques include 
different filtering approaches, such as the linear equalizer (a filter approximating 
the inverse of the charmel) and the decision feedback equalizer (DFE). These 
can be extended to the array signal case. Another often employed algorithm is 
the maximum likelihood sequence estimator (MLSE). The latter is often 




arcsin(a.) 
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regarded as being optimal, as it is derived from the maximum likelihood 
principle. 

As described above, the conventional DWILSP algorithm acts as a spatial 
diversity combiner, collecting the spatially spread energy in an efficient way. 
Thus, it would be desirable in the case of a channel spreadmg the transmitted 
energy both in space and time, to design an algorithm that performs diversity 
combining jointly over space and time. Such algorithms have been proposed but 
are unfortunately quite complex, both with regard to computational aspects and 
detection procedures involved. Also, these conventional approaches require an 
oversampling of the received signal. However, based on the well-known RAKE- 
approach, a space-time algorithm according to the present invention can be 
derived with the DWILSP algorithm as its elementary building block. 

Regarding the source signals as temporally white, the time dispersive case 
can be reformulated according to the frequency flat data model. The DWILSP 
algorithm can then be adopted to estimate different time-arrivals separately. This 
step thus performs spatial combing. Then, the different time-arrival estimates 
are combined temporally. This technique according to the present invemion thus 
constimtes a RAKE-combiner, exploiting both the spatial and temporal structure 
of the measured array signal, as well as the finite alphabet property of the 
modulated source signal. Moreover, this novel technique provides high 
performance at a low computational complexity, while at the same time lending 
itself to a simple and straightforward implementation. 

The approach taken here is based on estimation of different time arrivals 
of the desired user signal separately, instead of trying to invert or equalize the 
filter representing the channel. A final estimate is achieved by a combination of 
the estimates of the different time arrivals. 

To provide a foundation for understanding diversity combining according 
to the present invention, traditional RAKE techniques are fu-st described. The 
RAKE combiner was originally proposed for direct sequence spread spectrum 
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(DSSS) systems operating on time-dispersive channels. Consider the data model 
for a frequency selective channel: 



^(n) = ]CM«-*7;)+n(n) (22) 



Where s is a DSSS signal. A significant property of DSSS signals is that they are 
5 wideband signals. The wideband property is achieved by spreading the original 
data sequence with a high rate spreading code, whose elements are called chips, 
each with a duration of T, seconds. Each original data symbol thus contains 
several chips, and the spreadmg code is designed to have an autocorrelation 
function resembling white noise, such that symbols shifted more than one chip 

10 length apan are approximately uncorrelated. This type of signal is commonly 
used, for example, in radiocommunication systems that operate in accordance 
with code division multiple access (CDMA) techniques. • 

The DSSS RAKE combiner estimates each time-arrival s(n-kTc) by 
exploiting the autocorrelation property of the spreadmg sequence. The L-Hl 

15 signal estimates are then temporally combined to yield a finite signal estimate. 
The total scheme is thus equivalent to an L+1 order diversity combiner (if the 
channel taps h^ are uncorrelated). The conventional RAKE combiner can be 
illustrated as in FIG. 8, where each block 800 provides a time delay T, and the 
multiplication by c(n) at each multiplier 810 represents the despreading 

20 operation. The temporal branches seen in FIG. 8 are often referred to as 

"RAKE fingers" but are referred to herein as "spatio-temporal signal estimates" 
when used to refer to branches of a modified RAKE combiner wherein the 
DWILSP algorithm is used to provide for spatial combination. The ouq)uts of 
each RAKE finger are then temporally combined at block 820 by a diversity 

25 combining technique as will be described below. 
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Now consider the array signal model for the single user, frequency 
selective case: 



^("> = E V(«-^7;)+«(n) (23) 



The RAKE approach can also be applied to the array (unspread) signal 
case. Instead of spread symbols, as in the DSSS data model of equation (22), 
consider blocks of symbols. If the user signal is sufficiently temporally white', 
shifted versions, by an amount of seconds or more, become approximately ' 
uncorrelated. A block of symbols thus acts as the spreading sequence in the 
DSSS case, and different time-arrivals can be viewed as different user signals in 
the frequency flat case. Then the DWILSP type algorithm can be used to 
estimate the different time arrivals separately. 

The RAKE approach can be generalized to the multi-user case. 
Considering the different time arrivals as different signals, the double sum in the 
multi-user model can be rewritten according to equation (22), corresponding to 
15 the frequency flat case with d(L + 1) users: 



x(n) = i:V»^v(n) ^24) 



Thus, the RAKE combiner for the array signal case is shown in FIG. 9. where 
delayed versions of the received symbols are provided by blocks 900. The 
despreading operation in HG. 8 is replaced by the conventional DWILSP 
algorithm in blocks 910 which provide spatio-temporal signal estimates to the 
temporal combining block 920. 

Applying the conventional DWILSP algorithm to estimate time arrival k 
in the model of equation (23). it is seen that with a known temporally white user 
signal, the estimate of h, is still consistent. However, the filter tap correlations 



BNSOOCfD: <WO_9844655A2J_^ 



SUBSTITUTE SHEET (RULE 26) 



wo 98/44655 PCT/SE98/00S18 



-21- 

introduce an ISI tenn into the signal estimate in a manner similar to that 
described above. 

To overcome this problem associated with prior RAKE combiner efforts, 
i.e., to mitigate the effect of filter tap correlation, die self-synchronizing 
5 techniques described in die above exemplary embodiments can be applied to 
provide the spatio-temporal signal estimates as shown in FIG. 10. Therein, the 
delay blocks 1000, spatio-temporal signal estimators 1010 and temporal 
combining logic 1020 operate as described above. 

In the case of DSSS, maximum ratio combining (MRC) is often employed 
10 as the temporal combining technique employed in block 820 of FIG. 10. MRC 
maximizes die output signal-to-noise ratio (SNR), given independent noise in 
each finger and uncorrelated filter taps. The combined signal estimate is given 
as: 

(25) 

15 where Sk(n) is the output of the k:di finger, A/ die conjugate of die 
corresponding filter tap, and o^ is the finger noise variance. 



Modified Maximum Ratio Comhininy 

As described in the above incorporated by reference article entitled 
"Decoupled Separation of Digitally Modulated Signals Arriving at an Antenna 
Array", the conventional MRC approach can be modified with the conventional 
DWILSP algorithm. For example, before projection onto die finite symbol 
alphabet, die k:th estimate of die i:di user signal after the final iteration can be 
expressed as: 
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W = a/*ytt(«)+Ptt(«) (26) 

where represents a small bias (usually negligible), and 6^(n) is a noise 
term due to scaled thermal noise v(n) plus cochannel and self interferences s,(n). 
(j^i) OR (l^ic). This noise term can. with good accuracy, be considered L 
temporally white Gaussian and the noise in different signal estimates are 
approximately uncorrelated, i.e. E[Bj,fiJ=0. for Q^i) OR (l;^k). 

Ignoring the bias a, each signal estimate of Si(n) is automatically 
normalized in amplitude (PSK: | s(n) | =1) by DWILSP. Funhemiore, as the 
noise terms B(n) are uncorrelated between branches, one way of combining the 
estimates would be to use MRC at equation (25) as shown below: 



S,-(/i) = Proj 



(27) 



15 



20 



In equation (27), the operator (Proj) means projection onto the finite 
alphabet and o^^ is the variance of 13^(n). which can be estimated as Var(Proj(Sik)- 
§uc). For BPSK modulation. DWILSP projects symbols onto the alphabet of H-/-1 
and only the variance of the real part of the noise should be considered. 

However, when DWILSP fails completely at estimating a time-arrival of a 
signal, for example when a tap h^ in equation is weak or doesn't exist, the false 
signal estimate is stiU noimalized in amplitude, resulting in noise saturation 
effects, i.e. the variance a" of B(n) reaches a limit. To suppress bad 
estimates/fingers, simulations have indicated that a better weighting than ordinary 
MRC is: 
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^ 1 



(28) 



10 



Note that the value of the exponent (i.e., 4) in equation (28) is not 
critical. Any value in the range 3-6 results in approximately the same bit error 
rate (BER). 

Temporal Comhinino Using DWTLSP 

Although the DWILSP algorithm was originally intended for applications 
with antenna arrays, Applicants have recognized that tiiis technique is also a 
general diversity combiner. Consequently, the conventional DWILSP technique 
can also be used for the temporal combining function depicted in blocks 820, 920 
and 1020 of FIGS. 8-10, respectively. These exemplary embodiments of the 
present invention are illustrated in FIG. 11, wherein blocks 1100 provide delayed 
samples, blocks 1110, which are labelled simply as estimators to reflect a generic 
inclusion of traditional RAKE fingers, conventional DWILSP spatio-temporal 
estimators or usage of the self-synchronizing techniques according to the present 
invention, provide signal estimates and block 1120 shows the DWILSP technique 
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ftinctioning as the temporal combiner. Used in this way, stacking the L+1 (soft) 
estimates s^ of user signal i, a matrix equation is obtained as: 



= ws+Q (29) 
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15 



20 



where the left hand side of the matrix corresponds to the measured array 
signal X, s is the true signal and Q is a noise term. The column vector w can be 
interpreted as a temporal channel vector, representmg the delay profile of the 
channel. The column vector w also has a direct correspondence to the combining 

weights h^/al in MRC, or 1/a^ m the modified scheme, as solving 

equation (29) using DWILSP is essentially a search for the best diversity 
combining weights. Note that with the DWILSP algorithm employed for 
temporal combining, as well as to provide the spatio-temporal signal estimates, 
the finite alphabet property is used twice. 

There are several benefits associated with using the DWILSP technique 
for temporal combining in a RAKE receiver. For example, DWILSP is very 
robust in cases where not all diversity channels contain the signal of interest. 
Also, there is no specific assumption made regarding noise color, and the amount 
of fading correlation. Regardless of the temporal combining scheme applied, 
however, tracking of the combining weights for blocks/bursts of data, is 
desirable as it provides valuable information about time variations of the effective 
channel length. Also, this information is usefiil for frame synchronization in 
TDM A systems. 

To summarize the novel diversity combining techniques according to the 
present invention, an exemplary series of steps is illustrated by way of the 
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flowchart of Figure 12. First, at block 1200, the known training sequence (e. 
the CDVCC in D-AMPS) is used to obtam synchronization, and estimate the 



synchronizing technique described above can be used to obtain a signal estimate 
Sac of time path k at block 1210. Next, the signal estimates can be temporally 
combined at block 1220, by either (1) estimating the variance of 0-^(1) and using 
modified MRC according to equation (28) or using the conventional DWILSP to 
perform temporal combining. 

The performance of RAKE receivers according to the present invention 
was evaluated numerically at two different settings of a 5-tap FIR channel. The 
local scatter model was used to model each filter tap. The filter taps were 
modelled as statistically independent, therefore and the standard version of the 
DWILSP algorithm was used to provide the spatio-temporal signal estimates. 

In the simulations, the receiving antenna was chosen to be a 10-element 
uniform linear array (ULA). Three equally powered cochannel users were 
placed at nominal DOA:s [30^ 0^ -45^], relative to the array broadside. BPSK 
data was transmitted in bursts of 150 bits. Each burst included a 15-bit m- 
sequence, periodically extended to 19 bits, which was used as a training sequence 
for initialization of the receiver algorithm. 

As a comparison, the performance of the conventional MLSE technique 
was also evaluated. The MLSE was run twice: once with the exact chamiel and 
interference covariance matrix as a benchmark, and also using maximum 
likelihood estimates of these parameters obtained from the training sequence to 
provide a more realistic evaluation of MLSE performance. The MLSE was 
implemented with die Viterbi algorithm. 

In the first simulation, the relative average power in each tap was the 
same. Assuming a larger angular spread for the late arrivals (but the same 
nominal DO A), the cluster width standard deviations for the different taps were 
[2**, 3**, 4"*, 5*", e'']. The total signal power is here defmed as the sum of the 



filter length L. Then, initialized with the training sequence, the self- 
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powers in each filter tap. The resulting BER for the user at 0=0« as a function 
of the element Eb/N,, is shown in FIG. 13. 

Therein, staning from top right, the first curve shows the performance of 
the MLSE run with estimated channel parameters and interference covariance, 
5 and it is seen that this method has a performance that is limited by the cochannel 
interference. This characteristic is not seen in the other curves. The second 
curve shows the performance of the RAKE-combiner using standard MRC 
temporal combining. Moving to the third curve from the top right, it is seen that 
performance has been improved by about 2dB using modified MRC. Another 

10 2dB is gained by employing DWILSP combining. The last curve shows the 

performance of the MLSE run with the exact parameters. Considering the fact 
that the MLSE is run with estimated parameters in a practical application, e.g. 
the first curve, the RAKE-combiner using DWILSP as the temporal combining 
yields very good performance. 

15 In the second simulation, FIG. 14, the channel setting was adjusted to 

simulate a hilly terrain environment. The mean filter tap powers were set as [0, 
0, -20, -20, -6] dB and die cluster width standard deviations were [2°, 3", 6°, 
2°]. Each tap corresponded to a direction of [0", 1°, Q", 0°, io°] relative to the 
nominal DOA's, i.e., the late arrival impinging from a somewhat different 

20 direction compared to the early ones. The relative performance of the different 
algorithms resemble the results from FIG. 13. But this time, with an easier 
channel setting, the difference in performance between the DWTLSP-RAKE and 
MLSE using estimated parameters is larger. In fact, the DWILSP-RAKE comes 
close to the MLSE run with the exact parameters at high signal-to-noise ratios. 

is, of course, possible to embody the invention in specific forms other 
than those described above without departing from the spirit of the invention. 
The embodiments described above are merely illustrative and should not be 
considered restrictive in any way. The scope of the invention is determined by 
the following claims, rather than the preceding description, and all variations and 
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equivalents which fall within the scope of the claims are intended to be embraced 
therein. 
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WHAT IS CLAIMED IS: 

1- A method for processing a signal received in a 
radiocommunication system comprising the steps of: 

(a) sampling said received signal, whereby intersymbol 
5 interference can be introduced into said received signal; 

(b) generating at least one channel estimate associated with said 
received signal; 

(c) estimating said introduced intersymbol interference; 

(d) updating said at least one channel estimate associated with said 
10 received signal using said estimated introduced intersymbol interference; 

(e) repeating steps (b), (c) and (d) for a predetermmed number of 

iterations; and 

(f) outputting said received signal. 

2. The method of claim 1, wherein said step of sampling further 
15 comprises the step of: 

nonsynchronously sampling said received signal, 

3. The method of claim 1, wherein said step of repeating further 
comprises: 

performing, as said predetermined number of iterations, iterations 
20 until a value associated with said received signal converges. 

4. The method of claim 1 wherein said as least one channel estimate 
is generated using a decoupled weighted iterative least squares with projections 
(DWILSP) algorithm. 

5. The method of clakn 1, wherein said step of estimating further 
25 comprises the step of: 
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estimating said introduced intersymbol interference based on a type 
of modulation used to modulate said received signal. 

6. A system for processing a radio signal comprising: 
means for sampling said radio signal, whereby intersymbol 

5 interference can be introduced into said received signal; 

means for generating at least one channel estimate associated with 
said received signal; 

means for estimating said introduced intersymbol interference; 
means for updating said at least one channel estimate associated 
10 with said received signal using said estimated introduced intersymbol 
interference; 

means for iteratively applying said means for estimating, said 
means for generating and said means for updating for a predetermined number of 
iterations; and 

means for outputting said received signal. 

7. The system of claim 6, wherein said means for sampling further 
comprises: 

means for nonsynchronously sampling said radio signal. 

8. The method of claim 6, wherein said step of repeating- further 
20 comprises: 

performing, as said predetermined number of iterations, iterations 
until a value associated with said received signal converges. 

9. The method of claim 6 wherein said as least one chaimel estimate 
is generated using a decoupled weighted iterative least squares with projections 

25 (DWILSP) algorithm. 
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lo. The system of claim 6. wherein said means for estimating further 
comprises: 

means for estimating said introduced intersymbol interference 
based on a type of modulation used to modulate said received signal. 

11. A RAKE receiver comprising: 
an input node for receiving signal samples; 

a plurality of delay devices for generating delayed versions of said 
signal samples; 

a plurality of estimation branches, each receiving one of said 
delayed versions of said signal samples, and each for estimating a sample value 
therefrom; and 

a temporal combining unit for receiving said estimated sample 
values from said plurality of estimation branches, wherein said temporal 
combining unit uses a decoupled weighted iterative least squares widi projections 
15 (DWILSP) algorithm to combine said estimated sample values. 

12. The RAKE receiver of claim 11, wherein at least one of said 
plurality of estimation branches uses said DWILSP algorithm to estimate its 
respective sample value. 

13. The RAKE receiver of claim 11. wherem at least one of said 
20 plurality of estimation branches uses a modified version of said DWILSP 

algorithm which mitigates intersymbol interference introduced by filter taps to 
estimate its respective sample value. 

14. A RAKE receiver comprising: 

an input node for receiving signal samples; 

2^ a plurality of delay devices for generating delayed versions of said 

signal samples; 
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a plurality of estimation branches, each receiving one of said 
delayed versions of said signal samples, and each for estimating a sample value 
therefrom, at least one of said estimation branches using a modified version of a 
decoupled weighted iterative least squares with projections (DWILSP) algorithm 
5 to estimate its respective sample value, wherein said modified version operates to 
mitigate intersymbol interference introduced by filter taps; and 

a temporal combining unit for receiving and combining said 
estimated sample values from said plurality of estimation branches. 

15. The RAKE receiver of claim 14, wherein said temporal combining 
10 unit uses said DWILSP algorithm to combine said estimated sample values. 

16. The RAKE receiver of claim 14, wherein said temporal combining 
unit uses a maximal ratio combining technique to combine said estimated sample 
values. 
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